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Abstract 

We formulate and solve a model of two planes of antiferromagnetically cor- 
related electrons coupled together by a weak antiferromagnetic interaction of 
strength A. We show that in-plane antiferromagnetic correlations dramati- 
cally enhance the pairing effect of the interplane interaction. For the case 
where the in-plane correlation length find that the interac- 

tion A leads to spin pairing at a temperature T* ~ A, much higher than the 
usual BCS result exp{—J/X). We suggest that this is a possible explanation 
of the spin gap effects observed below T* ~ ISOi^T in Y Ba2Cu^0Q,Q. 
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Typeset using REVT^ 



It was recently argued that superconductivity and spin gaps in bilayer copper oxides such 
as Y Ba2Cu^0Q^x may be due to interplane pairing caused by the antiferromagnetic 
spin-spin interaction between the planes. Effects of this interaction have been observed in 
neutron scattering experiments on Y Ba2Cu^0Q,(i |Q. The high-Tc materials are also believed 
to have strong in-plane antiferromagnetic fluctuations. An alternative mechanism for spin 
gap formation in copper oxide materials, based on a single-plane theory of bosonic spin- 
waves, has also been discussed In this paper we determine the effect of the in-plane 
fluctuations on the interplane pairing interaction discussed previously. We find that they 
strongly enhance the between-planes interaction at wave vectors near the wave vector Q 
where the in-plane spin susceptibility is peaked. Taking into account this enhancement and 
the modification of the electron spectrum by the spin fluctuations we obtain an estimate 
for the onset temperature for the spin gap which is of the correct order of magnitude. 

Several different cases arise which we discuss in detail elsewhere 0. One issue is the 
relation between vector Q and the shape of the Fermi surface of the fermions: the vector 
Q might be a chord of the Fermi surface, its diameter and be larger than 2pp. [| In this 
paper we consider only the chord case. Another issue concerns the strength of the spin 
correlations. Here we assume that the spin system in each plane is very close to a T = 
critical point |§| resulting in long range antiferromagnetic fluctuations with a correlation 
length proportional to a power of temperature. A third issue is the nature of the fermionic 
excitations. One may distinguish the "spin liquid" case with spin-charge separation and 
fermionic spin excitations and the "Fermi liquid" case, where there is no spin-charge 
separation. Formally, the difference between these two pictures originates from the presence 
of an additional low energy mode ( gauge field ) in the spin liquid case , which results in 
a large relaxation rate for the fermions (so the electron propagator is (/5e^/^ — vf\p — Pf\)~^ 



all these conditions have trivial generalization for a non-spherical Fermi surface, for brevity we 
discuss only the circular case here. 
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PI). We treat both cases here. It is also beheved that in underdoped high-Tc materials 



electrons cannot tunnel coherently between planes |T2[, so we shall assume that all low 
energy excitations are confined to a plane. 

To model one plane of antiferromagnetically correlated fermions we write 

Hf = Y1 (^lAv)Cpa + Ji^i^-q (1) 

P Q 

where e(p) = t'^dp] — Pf) is the fermion dispersion near the Fermi surface, Sj = cl^Ca/sCpa^- 
It is also convenient to introduce the fermion mass m = pp/vp. The interaction Jq causes 
antiferromagnetic correlations peaked at wave vector Q. For definiteness we treat the inter- 
action in the RPA approximation and assume that the parameters are such that the spin 
susceptibility ^) is given by 

X{k,uj) = - r (2) 

/-I 

(3) 



+ (k-Q)2 + |a;|/r 

where xo is the susceptibility of the non-interacting fermions, k, the inverse correlation 
length, is assumed to be small and F is a 'microscopic' frequency scale. Presumably F ~ 1/m, 
or F ~ J/p|.. To fit the Cu NMR relaxation rates in high materials at high temperatures 
it is necessary to take = MT where M is a constant. We emphasize that although we 
have used the RPA to explain the form of @, this form is more general than the explanation 
and so are the following results which depend on @ only. The specific form of (|^) holds 
only if the wave vector Q < 2pp, so that at all wave vectors near Q a particle hole pair is 
available to damp the spin excitation. 

In the following we choose polar coordinates on the Fermi surface so that the points on 
the Fermi surface connected by Q correspond to angles ±6o- The form ([1|) applies to both 
the 'spin liquid' and Fermi liquid cases. In Q = 2pp case the functional form of x is different 
and depends on whether the fermion damping is small or large. 

We assume that the only coupling between different planes is an antiferromagnetic in- 
teraction between spins: 



i7.„, = AEsf^Sf\ (4) 

i 

where indices 1 and 2 distinguish planes in a bilayer and A is an interaction constant, 
assumed small. Neutron measurements imply that A ~ 200i^r, but certainly A -C J 
where J ~ ISOOfT is the exchange constant in one plane. 

The interaction leads to antiferromagnetic correlations between planes which we 
assume to be weaker than the in-plane correlations. An arbitrary weak A has also been 
shown to lead to a singlet pairing of spin excitations in different planes. In this work the 
antiferromagnetic correlations within each plane were not taken into and the temperature 
at which the spin pairing occurred was found to be very low (T^ ~ epe^^^'^^)- Here we show 
that in the presence of antiferromagnetic correlations the pairing interaction becomes much 
stronger at wave vectors near Q, the temperature at which the pairing occurs is greatly 
enhanced and the gap function becomes very anisotropic, opening first in a small region ( 
about {6 — Oq) ~ h/pf ) around the points connected by the vector Q and dropping rapidly 
away from these points as l/{6 — Oq)^. 

The physical argument is that because the susceptibility in one plane is very large at 
wave vectors near Q, a fermion at this wavevector polarizes the electrons in the neighboring 
plane in a large area around itself. Mathematically, we must construct the pairing vertex 
connecting a particle in one plane to a particle in the other. For small A this vertex will 
be linear in A and will be dressed by spin fluctuations in each plane: within RPA we have 
found that the dominant contribution to the dressed vertex V{k, uj) is that shown on fig 1 
which leads to 

V{k,u) = XJ'Qa-Wk,u;) 

where a is the lattice constant. Other contributions are negligible. To calculate the onset of 
the pairing from Eq. we must sum the ladder diagrams shown in Fig. 2. It is important 
to use the full Green's function, including the self energy due to spin exchange within one 
plane. This self energy has been studied by many authors. An approximation convenient 
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for our purposes is |T3|: 



V/ m aQ\uj\Jm 

J:{uj, 6) = , (6) 



where aq is a function of the order of unity if Q is away from 2pF but which diverges as 
Q ^.pp. We have verified that this formula apphes also in the spin-liquid case. 

It is important that the gap is due to the pairing of spin excitations on different planes, so 
the interaction (|^), although large in some region of momentum space, does not lead to self 
energy parts or vertex corrections. The gap equation, thus, follows from the summation of 
the ladder series in Fig. 2. By performing the ladder sum and integrating over the momenta 
in the direction normal to the Fermi surface we get 

^ [\uj\/T + piie^ + 6'^ + 2uee') + K^]^J[uj + i:{uj)]^ + A{e + uj,e'y 



where u = 1 — Q'^/{2ppY and we have set 9q = 0. The integration over the perpendicular 
momenta was possible because the main contribution to this integral comes from a narrow 
range near the Fermi surface {6p' ~ T/vp) where the interaction V{p — p', uj) does not vary 
significantly. 

To find the onset temperature we linearize (|^ and introduce scaled variables x and y via 
9 = kx/pf and 6' = ny/pp. The resulting equation is 



A ^ /■ , A,,+i{y)^l+y^ + ^\n + l + ^\ 



^"^"^^ 2aQMTa^J^J'^^\n + l + l\{l + y^ + x^ + 2uxy + ^\l\y 
where / and n are integers. From @ it is evident that A depends only on 1 + which 
enters only in the denominator of the kernel. Thus, A{6) is peaked about 6 = with a 
width K and decays for large 6 as 1/^^, and is similarly peaked about the lowest Matsubara 
frequency Un = ttT with the width F/?^ ~ T. The dimensionless kernel in (§) presumably 
has a largest eigenvalue w ~ 1, so T* is given by 

^ ^A 

2aQMa^J ^ ' 

Thus, apart from numerical factors the onset temperature T* is given by the bare interplane 
coupling constant A. For T <^ T* we may replace the sum over frequencies in by 



an integral; this integral is dominated by frequencies of the order of the zero temperature 
spin gap A(0) = A*; similarly, we must replace /t^ by MA* because low frequency spin 
correlations near the antiferromagnetic wavevector Q are eliminated by the spin gap. The 
result is that up to numerical factors T* in (|^) is replaced by A*. The gap takes its maximum 
value for angles 9 ^ V MA*; for larger 6 it given by A{6) ~ A* 



We emphasize that due to a strongly peaked and temperature dependent effective in- 
teraction, the pairing temperature and the gap scale as the interaction constant, unlike the 
usual BCS case where they are exponentially small. Although we have assumed specific 
form for the spin susceptibility @ with a temperature dependent correlation length, this 
assumption is not essential to our results. The enhancement of the between-planes pairing is 
due to the strong temperature dependence of J2qX'{lj^ = 0)^. This quantity is measurable 
via NMR T2 experiments and has been found to be large and strongly temperature 
dependent in FSaaCwsOr-^ [pHlSj. 

The development in the spin liquid case is essentially identical. The momentum inte- 
grated Green's function square is [/?|co'p/'^ — S(ct;, 9)]^^ but /^Icup/^ is still negligible compared 
to T,{uj,9) for the frequencies and angles of interest, so Eqs are not changed. 

We now consider the experimental implications. The pairing mechanism is much weaker 
in La2-xSrxCu04 because the antiferromagnetic interaction between Cu ions in different 
planes is frustrated, so that in tetragonal crystals Eq (H) becomes 

^f„* = AEsr^Sa (10) 

where 6 labels the four Cu sites in plane 2 equidistant from site i of plane 1. Eq (|10|) implies 
that V{k,u) in eq. (^) becomes 

V^^'ik, 00) = V{k, uj) cos{kj2)cos{ky/2) (11) 

Thus the singularity in the interaction is eliminated for commensurate spin fluctuations 
{krc, ky ~ vr) in tetragonal crystals. For orthorhombic crystals or for incommensurate spin 
fluctuations the singular part of the interaction is of order the square of the orthorhombicity 
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or incommensurability, and is therefore small. This is consistent with the observation that 
the spin gap opens at much lower temperatures in La2-xSrxCuOii than in Y Ba2Cu^0^j^x- 
In a Fermi liquid system with no spin-charge separation the opening of the spin gap 
implies that the material has become superconducting. In a spin liquid system, true su- 
perconductivity will only occur at a lower temperature where the charge carriers bose con- 
dense. The former scenario is consistent with the behavior of optimally doped or overdoped 
Y Ba2CuzO^+x and with La2-xSrxCu04 at all dopings, while the latter scenario is consis- 
tent with the behavior of underdoped YBa2Cu30Q+x- For example, in YBa2Cu30Qs, spin 
gap effects are observed in NMR below T* ~ 1507^ while the superconducting Tc ~ 60K. 
As previously pointed out [|[] there is also optical evidence |]16| for the existence of a gap 
above T^. The small value of the specific heat jump at in YBa2Cu30Qs 
sistent with this scenario. However, none of these observations (except the qualitative one 
that the spin gap opens significantly above Tc only in underdoped bilayer materials such as 
YBa2CusOQs) distinguishes the mechanism we have proposed from other possible origins 
of the spin gap. 

There is one qualitative disagreement with experiment. Because the gap opens first and 
is largest at the points on the Fermi surface connected by the wavevector where x{k,uj) is 
peaked, the low frequency antiferromagnetic spin fluctuations are suppressed more strongly 
than spin fluctuations at other wavevectors. In the high-T^ materials it is believed that 
the antiferromagnetic fluctuations are responsible for the enhancement of the Cu relaxation 
rate over the relaxation rates of the other nuclei therefore in our scenario the copper 
relaxation rate would drop more rapidly than the oxygen or yttrium rates as the spin gap 



opened, in apparent disagreement with experimental data on YBa2Cu30Qs [18 



Note added: As this manuscript was being prepared we learned that M. Ubbens and P. 
A. Lee [p9| have obtained results very similar to ours. 
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FIGURES 

FIG. 1. Diagrammatic representation of dominant contribution to pairing interaction. Here 
the light dashed line represents the intcrplanc interaction and the wavy line represents the dressed 
spin-fluctuation interaction between electrons in one plane. 

FIG. 2. Ladder sum leading to gap equation. Here the shaded rectangle is the interaction V 
defined in Fig 1. Note that electron lines are dressed by in-plane fiuctuations. 
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